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1. INTRODUCTION 

Dusty tori around active galactic nuclei (AGNs) play 
an important role in the classification of Seyfert galaxies. 
(Antonucci & Miller 1985; see also Antonucci 1993 for a re- 
view) . Seyfert galaxies observed from a face-on view of the 
torus are recognized as type 1 Seyferts (Sis) while those 
observed from an edge-on view are recognized as type 2 
Seyferts (S2s). Therefore, physical properties of dusty tori 
are of great interest. We briefly introduce three statistical 
studies investigating properties of dusty tori; 1) physical 
sizes of dusty tori based on water-vapor maser emission 
(Taniguchi & Murayama 1998), 2) ionization condition of 
the inner wall of tori based on high-ionization emission 
lines (Murayama & Taniguchi 1998a, b), and 3) viewing 
angle toward dusty tori based on mid-infrared color (Mu- 
rayama, Mouri, & Taniguchi 2000). Please see references 
for detailed discussion. 

2. DUSTY TORI OF SEYFERT NUCLEI POSED BY THE 
WATER VAPOR MASER EMISSION 

2.1. Water Vapor Maser Emission in Active Galactic 

Nuclei 

The recent VLBI/VLBA measurements of the H2O 
maser emission of the nearby AGNs, NGC 1068 (Gal- 
limore et a!. 1996; Greenhill et al. 1996; Greenhill & 
Gwinn 1997), NGC 4258 (Miyoshi et al. 1995; Greenhill 
et al. 1995a, 1995b), and NGC 4945 (Greenhill, Moran, 
& Herrnstein 1997), have shown that the masing clouds 
are located at distances of ~ 0.1 - 1 pc from the nuclei. 
These distances are almost comparable to those of molec- 
ular/dusty tori which are the most important ingredient 
to explain the observed diversity of AGN (Antonucci & 
Miller 1985; Antonucci 1993). It is therefore suggested 
that the masing clouds reside in the tori themselves (e.g., 
Greenhill et al. 1996). Therefore, the H2O maser emis- 
sion provides a useful tool to study physical properties of 
dusty tori which are presumed to be the fueling agent onto 
the supermassive black hole (cf. Krolik & Begelman 1988; 
Murayama & Taniguchi 1997). 

2.2. A Statistical Size of the Dusty Tori Inferred from 

the Frequency of Occurrence of H% O Masers 

The recent comprehensive survey of the H2O maser 
emission for ~ 350 AGNs by Braatz et al. (1997; hereafter 
BWH97) has shown that the H2O maser emission has not 
yet been observed in Sis and that the S2s with the H2O 
maser emission have the higher H I column densities to- 
ward the central engine. It is hence suggested strongly 
that the maser emission can be detected only when the 
dusty torus is viewed from almost edge-on views. This is 
advocated by the ubiquitous presence of so-called the main 
maser component whose velocity is close to the systemic 



one whenever the maser emission is observed because this 
component arises from dense molecular gas clouds along 
the line of sight between the background amplifier (the 
central engine) and us (see, e.g., Miyoshi et al. 1995; 
Greenhill et al. 1995b). 

Since the high H I column density is achieved only 
when we see the torus within the aspect angle, <f> = 
tan _1 (/i/26) (see Figure 1), we are able to estimate b 
because the detection rate of H 2 maser, P mase r, emis- 
sion can be related to the aspect angle as, -P maS cr = 
A ma scr/ (^mascr Y non _ maser ) = cos(90°-</>) where N maser 
and A non _ mascr are the numbers of AGN with the H 2 
maser emission and without the H2O maser emission, 
respectively. This relation gives the outer radius, b = 
h [2tan(90° — cos -1 P ro aser)] _1 - Table 1 shows that a typ- 
ical detection rate is P maS or — 0.05. However, this value 
should be regarded as a lower limit because some special 
properties of may be necessary to cause the maser emis- 
sion (Wilson 1998). If we take account of new detections 
of H2O maser emission from NGC 5793 (Hagiwara et al. 
1997) and NGC 3735 (Greenhill et al. 1997b) which were 
discovered by two other maser surveys independent from 
BWH97, the detection rate may be as high as ~ 0.1 (Wil- 
son 1998). Therefore, we estimate b values for the two 
cases; 1) P maS cr = 0.05, and P maS cr = 0.1. These two rates 
correspond to the aspect angles, ~ 2?9 and (j> — 5?7, 
respectively. In Table 2, we give the estimates of b for 
three cases, a — 0.1, 0.5, and 1 pc. If a > 1 pc, the 
H I column density becomes lower than 10 23 cm -1 given 
Mgas — 10 5 M Q . Therefore, it is suggested that the inner 
radius may be in a range between 0.1 pc and 0.5 pc for 
typical Seyfert nuclei. The inner radii of the H2O masing 
regions in NGC 1068, NGC 4258, and NGC 4945 are in- 
deed in this range (Greenhill et al. 1996; Miyoshi et al. 
1997; Greenhill et al. 1997a). We thus obtain possible 
sizes of the dusty tori; (a, b, h) — (0.1 - 0.5 pc, 1.67 - 8.35 
pc, 0.33 - 1.67 pc) for (j) ~ 5?7, and {a,b,h) = (0.1 - 0.5 
pc, 3.29 - 16.5 pc, 0.33 - 1.67 pc) for <f> ~ 2?9. All the 
cases can achieve A/hi > 10 23 cm -1 , being consistent with 
the observations (BWH97). 

3. HIGH-IONIZATION NUCLEAR EMISSION-LINE REGIONS 
ON THE INNER SURFACE OF DUSTY TORI 

3.1. High-ionization Emission Lines in Seyfert Galaxies 

Optical spectra of active galactic nuclei (AGN) show of- 
ten very high ionization emission lines such as [Fe VII], 
[Fe X], and [Fe XIV] (the so-called coronal lines). Accord- 
ing to the current unified models (Antonucci & Miller 1985; 
Antonucci 1993), it is generally believed that a dusty torus 
surrounds both the central engine and the BLR. Since the 
inner wall of the torus is exposed to intense radiation from 
the central engine, it is naturally expected that the wall 
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Fig. 1. — The geometry of the dusty torus. The torus is a cylinder of dust with a uniform density, characterized by the inner radius (a), 
the outer radius (b), and the full height (h). A half-opening angle of the torus is thus given as (f> = ta.n~ 1 (h / 2b) . 



Table 1 

a summary of the detection rates of the h 2 maser in active galactic nuclei studied by braatz, 

Wilson, & Henkel (1997) for the various samples 



Sample 


N 

1 v mascr 


^Vtotal 


-^mascr (%) 


Distance- limited 


All (S1+S2+L) 


15 


278 


5.4 


Seyfert (S1+S2) 


10 


198 


5.1 


S2 


10 


141 


7.1 


Magnitude- limited 


All (S1+S2+L) 


13 


241 


5.4 


Seyfert (S1+S2) 


8 


166 


4.8 


S2 


8 


112 


7.1 



Table 2 

Geometrical properties of the dusty tori inferred from the statistics of the H 2 maser emission 











-^mascr — 0.05 




P — 1 

1 mascr u - r 










<j> = 2?9 




= 5?7 


a (pc) 


h (pc) 


r h ot a (pc) 


b (pc) 


N m (cm" 2 ) 


b (pc) 


N m (cm" 2 ) 


0.1 


0.33 


0.43 


3.29 


3.3 x 10 24 


1.67 


6.5 x 10 24 


0.5 


1.67 


2.17 


16.5 


1.3 x 10 23 


8.35 


2.6 x 10 23 


1 


3.30 


4.30 


32.9 


3.3 x 10 22 


16.7 


6.5 x 10 22 



a The radius of the hot part in the torus; r hot = a + h. 
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can be one of the important sites for the HINER (Pier 
& Voit 1995). If the inner wall is an important site of 
HINER, it should be expected that the Sis would tend to 
have more intense HINER emission because the inner wall 
would be obscured by the torus itself in S2s. 

In order to examine whether or not the Sis tend to have 
the excess HINER emission, we study the frequency distri- 
butions of the [Fe VII] A6087/[O III] A5007 intensity ratio 
between Sis and S2s. The data were compiled from the 
literature (Osterbrock 1977, 1985; Koski 1978; Osterbrock 
& Pogge 1985; Shuder & Osterbrock 1981) and our own 
optical spectroscopic data of one SI (NGC 4051) and four 
S2s (NGC 591, NGC 5695, NGC 5929, and NGC 5033). In 
total, our sample contains 18 Sis and 17 S2s. The result 
is shown in Figure 2. It is shown that the Sis are strong 
[Fe VII] emitters than the S2s. In order to verify that this 
difference is really due to the excess [Fe VII] emission, we 
compare the [O III] luminosity between the Sis and S2s 
and find that the [O III] luminosity distribution is nearly 
the same between the Sis and the S2s (Figure 3). There- 
fore, we conclude that the higher [Fe VII] /[O III] intensity 
ratio in the Sis is indeed due to the excess [Fe VII] emis- 
sion rather than the weaker [O III] emission in the Sis. 
The presence of an excess [Fe VII] emission in Sis can 
only be explained if there is a fraction of the inner HINER 
that cannot be seen in the S2s. The height of the inner 
wall is of order 1 pc (Gallimore et al. 1997; Pier & Krolik 
1992, 1993). Therefore, given that the torus obscures this 
HINER from our line of sight, the effective height of the 
torus should be significantly higher than 1 pc. 

3.2. Three- Component HINER 

Although our new finding suggests strongly that part of 
the HINER emission arises from the inner walls of dusty 
tori, it is remembered that a number of S2s have the 
HINER. In fact, the fraction of Seyfert nuclei with the 
HINER is nearly the same between the Sis and the S2s 
(Osterbrock 1977; Koski 1978). If the HINER was mostly 
concentrated in the inner 1 pc region, we would observe the 
HINER only in the Sis. Therefore the presence of HINER 
in the S2s implies that there is another HINER component 
which has no viewing-angle dependence. A typical dimen- 
sion of such a component is of order 100 pc like that of the 
NLR. In addition, it is also known that some Seyfert nuclei 
have an extended HINER whose size amounts up to ~ 1 
kpc (Golev et al. 1994; Murayama, Taniguchi, & Iwasawa 
1998). The presence of such extended HINERs is usually 
explained as the result of very low-density conditions in 
the interstellar medium (nn ~ 1 cm~ 3 ) makes it possible 
to achieve higher ionization conditions (Korista & Ferland 
1989). 

The arguments described here suggest strongly that 
there are three kinds of HINER; 1) the torus HINER (r < 1 
pc), 2) the HINER associated with the NLR (10 < r < 100 
pc), and 3) the very extended HINER (r ~ 1 kpc). A 
schematic illustration of the HINER is shown in Figure 4. 

3.3. Dual- Component Photoionization Calculations for 

HINER 

Any single-component photoionization models under- 
prcdict higher ionization emission lines (see Murayama 
& Taniguchi 1998b and references therein). We therefore 
proceed to construct dual-component models in which the 



inner surface of a torus is introduced as a new ionized-gas 
component in addition to the traditional NLR component 
with the photoionization code CLOUDY (Ferland 1996). 
The single-cloud model suggests that the ionization pa- 
rameter lies in the range of \ogU ~ —1.5 2. As for the 

electron density, it is often considered that the inner edges 
of tori have higher electron densities, e.g., n c ~ 10 78 cm~ 3 
(Pier & Voit 1995). Because the largest [Fe VII] /[O III] 
ratio of the observed data is ~ 0.5, [Fe VII] /[O III] of the 
torus component must be greater than 0.5. However, we 
find that ionization-bounded models cannot explain the 
observed large [Fe VII] /[O III] values by simply increas- 
ing electron densities up to 10 9 cm~ 3 . Further, such very 
high-density models yield unusually strong [O I] emission 
with respect to [O III]. We therefore assume "truncated" 
clouds with both large [Fe VII] /[O III] ratios and little 
low-ionization lines for the HINER torus. The calcula- 
tions were stopped at a hydrogen column density when 
[Fe VII] /[O III] = 1. We performed photoionization calcu- 
lations described above and we finally adopted the model 
with jih = 10 7,5 cm~ 3 and logC/ = —2.0 representative 
model for the HINER torus with taking [Fe X]/[Fe VII] 
ratios predicted by the calculations into account. 

Now we can construct dual-component models combin- 
ing this torus component model with the NLR models. 
In Figure 5, we present the results of the dual-component 
models. Here the lowest dashed line shows the results of 
the NLR component models with a = — 1, logU = —2, as 
a function of nn from 1 cm~ 3 to 10 6 cm~ 3 . If we allow 
the contribution from the torus component to reach up to 
~ 50 % in the Seyferts with very high [Fe VII] /[O III] ra- 
tios, we can explain all the data points without invoking 
the unusual iron overabundance. Note that the majority 
of objects can be explained by simply introducing a ~ 10 
% contribution from the HINER torus. 

4. NEW MID-INFRARED DIAGNOSTIC OF THE DUSTY 
TORUS MODEL FOR SEYFERT NUCLEI 

4.1. The New MIR Diagnostic 

The current unified model of active galactic nuclei 
(AGNs) has introduced the dusty torus around the cen- 
tral engine (Antonucci 1993). Therefore, it is urgent to 
study the basic properties of dusty tori (e.g., Pier & Kro- 
lik 1992). Utilizing the anisotropic property of dusty torus 
emission, we propose a new MIR diagnostic to estimate a 
critical viewing angle of the dusty torus between type 1 
and 2 AGNs. 

Because of the anisotropic properties of the dusty torus 
emission, the emission at A < 10 /im is systematically 
stronger in type 1 AGNs than in type 2s while that at 
A > 20 ^m is not significantly different between type 1 
and type 2 AGNs. Therefore the luminosity ratio between 
3.5 /zm and 25 /nm is expected to be highly useful to dis- 
tinguish between type 1 and 2 AGNs (Figure 6). Here we 
define the above ratio as 

R = l0gJ/ 3 .5 Mm /j, 3 . 5Mm /^25Aim /i/ 25M m ■ 

4.2. Results & Discussion 

We adopt three samples chosen by different selection cri- 
teria and compiled photometric data in L, N, and IRAS 
25 fim bands: 
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Fig. 4. — A schematic illustration of the three-component model for the HIENR. Note that the torus HINER consists of many small ionized 
gas clumps like the clumpy HINER in the NLR. 




log([OI]6300/[OIII]5007) 



Fig. 5. — Dual-component photoionization models arc shown in the diagram of [Fe VII] /[O III] vs. [O I]/[0 III]. The circles are Sis, the 
triangles are SI. 5s, and the squares are S2s. The filled symbols denote the objects with [Fe X] emission, while the open symbols denote the 
objects without [Fe X]. The numbers labeling the lowest dashed line represent log ran- The percentages represent contribution of the HINER 
[O III] flux to the total [O III] flux. 
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Fig. 6. — Basic concept of our MIR diagnostic. The upper panel shows typical spectra of the torus emission for Sis (upper) and for S2s 
(lower). The lower panel shows how the 3.5 fim to 25 /im flux ratio yields the viewing angle toward the torus. 



1. 18 Sis and 6 S2s from the CfA Seyfert galaxies 
(Huchra & Burg 1992) 

2. 20 Sis and 4 S2s from the sample of Ward et al. 
(1987), which is limited by the hard X-ray flux from 
2 to 10 keV 

3. 11 Sis and 11 S2s from the sample of Roche et al. 
(1991), which is composed of TV-band bright objects 

Since some objects are included in more than one sample, 
there are 31 Sis and 14 S2s in total. 

The type 1 Seyferts are clearly distinguished from the 
type 2s with a critical value R ~ —0.6; R > —0.6 for 
type Is while R < —0.6 for type 2s (Figures 7a-d). If we 
apply the Kolmogrov-Smirnov (KS) test, the probability 
that the observed distributions of Sis and S2s originate in 
the same underlying population turns out to be 0.275 %. 

The upper panel of Figure 8 shows the theoretical mod- 
els of Pier & Krolik (1992, 1993), which are characterized 



by a (the inner radius of the torus), h (the full height 
of the torus), r r (the radial Thomson optical depth), r z 
(the vertical Thomson optical depth), and T (the effective 
temperature of the torus) [see Figure 9] . The intersection 
of each model locus with R = —0.6 gives a critical view- 
ing angle. The critical viewing angle is expected to be 
nearly the same as the typical semi-opening angle of the 
ionization cones observed in Seyfert nuclei, ~ 30°- 40° (cf. 
Lawrence 1991 and references therein). Figure 9 shows 
that only two models give reasonable critical viewing an- 
gles, ~ 46°- 50°though these values are slightly larger than 
the semi-opening angle of the cone. The model with a/h = 
0.3 may be suitable for tori in Seyfert nuclei because this 
inner aspect ratio gives a semi-opening angle of the torus, 
~ 30° , being consistent with those of the observed ionized 
cones. Although there are some contaminations from the 
host galaxies, circumnuclear starbursts, and dust emission 
in the narrow-line regions, the new diagnostic provides a 
powerful tool to study the critical viewing angle. 
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Fig. 7. — Histogram of R(L, 25) for the CfA Seyferts (a), the sample of Ward et al. (1987) (6), the sample of Roche et al. (1991) (c), and 
the total sample (d). Galaxies shown by white bars are likely to suffer from contamination and are not used in our analysis. 
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Fig. 8. — Upper panel: relationships between R(L, 25) and the viewing angle for six dusty torus models given in Table 4. Lower panel: 
distributions of the observed R(L, 25) values. Galaxies shown by open triangles are likely to suffer from contamination and are not used in 
our analysis. 
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Fig. 9.— 



Schematic illustration of the geometrical configuration of the dusty torus model. 



